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We develop a method of the ab-initio calculation for the core-level x-ray photoemission spec- 
troscopy (XPS). By calculating one-electron states in the presence of core-hole potential, excited 
states are constructed by distributing electrons on these one-electron states. The overlap integrals 
between the excited states and the ground state are evaluated by using the wavefunctions with and 
without the core-hole potential, and finally the XPS spectra are obtained. Implementing the proce- 
dure, we obtain the spin-resolved 3s-core XPS spectra in ferromagnetic iron without any adjustable 
parameters, in good agreement with the experiment. The spectral shapes are quite different for 
different spin channels. We explain the behavior in terms of the difference in the one-electron states 
screening the core-hole potential. 

PACS numbers: 79.60.-i 71.15.Qe 71. 20. Be 



INTRODUCTION 



Electronic structure calculations based on the den- 
sity functional theory with the local approximation 
have provided valuable information on the ground state 
properties.— Various attempts beyond the local density 
approximation (LDA) have also been worked out. These 
approaches have made great success on predicting ac- 
curate electronic structures without introducing any ad- 
justable parameters.— However, the LDA and its exten- 
sions have not been suited to treat highly-excited states, 
and are usually limited to calculate the ground-state 
properties. For example, although the x-ray photoemis- 
sion spectroscopy (XPS) from inner-core states is a pow- 
erful method to investigate the dynamical response to 
the local perturbation associated with the photo-excited 
core-hole,— the corresponding spectra have not been cal- 
culated from the first principles. 

The purpose of this paper is to develop a method of the 
ab-initio calculation for the inner-core level XPS spec- 
tra. It is known that the spectra show asymmetric shape 
as a function of binding energy near the threshold in 
metals A This phenomenon is related to the singularity 
in the response function near the Fermi edge^^ The 
spectra sometimes show extra peaks called satellites in 
systems of strongly-correlated electrons. One example 
is in systems of rare-earth-metal compounds, where the 
main and satellite peaks are explained as "well-screened" 
and "poorly-screened" states, respectively.-'- Another no- 
table example is a satellite around 6 eV separate from 
the threshold on the 2p-XPS in a ferromagnetic metal 
NiiiS It arises from a combined effect of screening the 
core-hole potential and the interaction between electrons. 
Feldkamp and Davis (FD) analyzed the 2p-XPS spectra 
in Ni by developing a numerical method on the Hub- 
bard model.™ Evaluating the overlap between the excited 
states and the ground state, they obtained the spectra in 
good agreement with the experiment. 



The method we develop in this paper is as an exten- 
sion of the FD theory to an ab-initio level. First, we 
carry out the band calculation within the LDA to ob- 
tain the one-electron states in the ground state. We use 
the augmented plane wave (APW) method instead of the 
Korringa-Kohn-Rostoker (KKR) method, since we need 
the wavefunctions later. Next, we consider a system of 
supercells with one core-hole per cell. The core hole is 
introduced by removing an electron from the core level, 
which wavefunction is assumed to be localized. In this 
system, we carry out the band calculation to obtain the 
one-electron states. To guarantee the charge neutrality, 
we add one extra electron in each supercell. To be ex- 
act, there should exist only one core-hole in crystal for 
the XPS event. In this situation, as known in the im- 
purity problem, the local charge neutrality is naturally 
satisfied according to the Friedel sum ruleji 2 - The super- 
cell method works better as increasing the cell size. Note 
that not only the effect of core-hole potential but also 
that of electron-electron interactions are taken into ac- 
count in the one-electron states thus obtained through 
the exchange-correlation potential. Discretizing the mo- 
mentum space into finite number of points, we distribute 
electrons into the one-electron states to construct the ex- 
cited states. The lowest-energy state with core hole is 
given by piling electrons into all one-electron levels up 
to the Fermi level. One electron-hole pair, two electron- 
hole pairs, and so on, are created from this core-hole 
state. Finally, we calculate the XPS spectra by evalu- 
ating the overlaps between these excited states and the 
ground states with the help of the one-electron wavefunc- 
tions. 

We substantiate this scenario by applying the pro- 
cedure to the 3s-core XPS in the ferromagnetic Fe 
metal . 13 ' 14 ] 15 ! 16 ! 17 ! 18 ! 19 The spectra are resolved with re- 
spect to the spin of photoelectron or equivalently the spin 
of core hole. Hereafter we call the majority and minor- 
ity spin as the up and down spin, respectively. The core 
hole with up (down) spin indicates that one core elec- 
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tron with up (down) spin is missing. The spectra show a 
strong dependence on the spin^ for the core hole with 
up spin, the spectral intensity is distributed in a wide 
region of binding energy with a clear satellite peak; for 
the core hole with down spin, the spectra have a slightly 
asymmetric peak with no satellite peak. These features 
indicate that the screening effect on the core hole with 
up spin is quite different from that with down spin. So 
far, such spectra have been analyzed using a single band 
Hubbard model,— but it is too simple to compare quan- 
titatively the experimental data.— We demonstrate that 
the present scheme works well by comparing the calcu- 
lated spectra with the experimental data. 

Since the core hole has a spin-degrees of freedom, the 
exchange effect plays an important role, which is different 
from the impurity potential problem in dilute alloys . 21 ' 22 
We could explain the strong spin dependence of the spec- 
tra in relation to 3d-bands changed by the core hole. For 
the core hole with up spin, the 3d electrons with down 
spin are strongly attracted by the core-hole potential, re- 
sulting in an overscreening, while the 3d electrons with 
up spin are little influenced by the core-hole potential be- 
cause of the exchange effect, resulting in an antiscreen- 
ing due to the repulsion between 3d electrons. The local 
spin moment is reversed at the core-hole site. On the 
other hand, for the core hole with down spin, the 3d elec- 
trons with up spin are strongly attracted by the core-hole 
potential, but the unoccupied 3d states available to the 
screening are limited, resulting in an incomplete screen- 
ing. The shortage of the screening is compensated by the 
3d electrons with down spin. The spin-dependent XPS 
spectra are interpreted in terms of these changes in the 
3d one-electron states. 

Note that the 3s satellite is conventionally interpreted 
as a multiplet splitting between the 3s and 3d electrons 
on the basis of the atomic model.— However, the elec- 
tronic structure given in this paper is quite different from 
what the atomic model predicts. The number of 3d elec- 
trons is increased nearly one at the core-hole site in accor- 
dance with Friedel sum rule, and the one-electron states 
are strongly modified from the ground state, and even 
quasi-bound states appear. 

Furthermore, we calculate the ls-core XPS spectra to 
clarify the situation by comparison with the 3s-core spec- 
tra. Since the Is wavefunction is so localized that the 
exchange effect on the 3d one-electron states is expected 
to be relatively small. This is actually the case, and we 
obtain the XPS spectra, which are less dependent on the 
spin. 

The present paper is organized as follows. In Sec. 2, 
we formulate the XPS spectra with the ab-initio method. 
In Sec. 3, we present the calculated XPS spectra for Fe 
3s- and ls-core in comparison with the experiment. The 
last section is devoted to the concluding remarks. 



II. FORMULATION 

A. XPS spectra 

We consider the situation that a core electron is excited 
to a high energy state with energy e by absorbing the x 
ray with energy ui q . Then the probability of finding a 
photoelectron with energy e may be proportional to 



2ttM 2 ^ \(f\s a \g)\ 2 S(oj q + E g -e- E f ), (1) 
/ 

where w represents the transition matrix element from 
the core state localized at a particular site to the state 
of photoelectron. The s CT is the annihilation operator of 
a relevant core electron, which is assumed to have only 
spin a as the internal degrees of freedom. State \g) rep- 
resents the ground state with energy E g , while |/) rep- 
resents the final state with energy Ef. For the defini- 
tion of |/), the photoelectron is excluded. In the actual 
calculation, we replace the (^-function by the Lorentzian 
function with the full width of half maximum (FWHM) 
2r s to taking account of the life-time broadening of the 
core level. It should be noted here that we can specify 
the spin of core hole by resolving the spin of photoelec- 
tron. Since the core-hole potential induces the screening 
in the final state, the electronic structure would become 
different from the ground state. Therefore, the spectra 
would be spin-dependent in the magnetic system, due to 
the exchange interaction. 



B. Overlap integrals 

The KKR method could solve a single core-hole prob- 
lem in the final state. Actually it is successfully applied to 
alloy problems with combined to the coherent-potential 
approximation (CPA) i 23 i 24 However, it may not be useful 
in the present problem, because we need the information 
of wavefunctions. Instead, we consider a periodic array of 
supercells containing one core-hole per cell, and calculate 
the one-electron states by means of the augmented plane 
wave (APW) method. The larger the unit cell size is, 
the better results are expected to come out. Inner-core 
states, for example Is, 2s, 2p, and 3s in Fe, are treated 
as localized states within a muffin-tin sphere, so that 
we could directly specify the core hole site. To ensure 
the charge neutrality, we assume n e + 1 band electrons 
per unit cell instead of n e band electrons in the ground 
state. One additional electron per unit cell would not 
cause large errors in evaluating one-electron states in the 
limit of large unit cell size. We write the resulting single 
particle states with energy eigenvalue e„(k) as 

ip n k{r) = — L= ^2 <?W(r - Rj) exp(ikRj), (2) 
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Figure 1: Sketch of a supercell containing core holes in bcc 
Fe. 

with 

0„ k (r-R J ) = U „ k (r)e*( r - R ^, (3) 



Table I: Screening charge with respect to the d symmetry 
within the muffin-tin sphere on the 3s-core-hole site. 

An d | An d | An d (total) 
up spin hole -1.44 2.38 0.94 
down spin hole 0.47 0.48 (L95_ 

III. CALCULATED RESULTS 

We apply the above method to calculate the XPS spec- 
tra for the 3s and Is cores in Fe metal. We use a supercell 
containing 27 iron atoms, as illustrated in Fig. [TJ The 
first Brillouin zone (BZ) for the supercells becomes much 
smaller than the original first BZ. Accordingly, the energy 
bands in the absence of core hole is folded. To demon- 
strate the reliability of our supercell calculation, we show 
in Fig.[2]the DOS calculated in the supercell system with- 
out core hole. The DOS resolved into the majority (up) 
spin and minority (down) spin reproduces well the previ- 
ous results.— In addition, the DOS is projected onto the 
d symmetry (g?-DOS) within the muffin-tin sphere, which 
curve is close to the total DOS. This indicates that the 
most part of the DOS arises from the 3d states. 



where u„k( r — Rj) is a periodic function with changing 
unit cells, and j runs over N c unit cells. We use these one- 
electron states as substitutes for the states under a single 
core-hole. We distribute N e (= n e N c ) band electrons 
on these states to construct excited states. In addition, 
we carry out the band calculation in the absence of the 
core hole with assuming n e band electrons per unit cell. 
The wavefunction and eigenenergy are denoted as i/^k ( r ) 
and e^^k), respectively. All the lowest N e levels are 
occupied in the ground state. Then, the matrix elements 
connecting the ground and final states are expressed as 



with 
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Figure 2: DOS calculated in a system of supercells with no 
core hole in ferromagnetic iron. The solid line represents the 
DOS projected onto the d-symmetry within the muffin-tin 
sphere. 



A. 3s-core hole 



where the integral is carried out within a unit cell. Sub- 
scripts i — (n, k) and i' — (n',k') are running over oc- 
cupied states in the final state and in the ground state, 
respectively. The corresponding energies are 



E f = ei 

TP - ^° 
- e i 
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(6) 
(7) 



Substituting Eqs. gj, ©, and into Eq. (JTJ) , we obtain 
the XPS spectra. 



1. one-electron states 

How are the one-electron states modified by the core- 
hole potential? We explain the change through the 
change in the g?-DOS at the core-hole site. Figure [3] 
shows the d-DOS calculated in the presence of the 3s-core 
hole, and table [J lists the screening charge Arid, which is 
defined by the change in the occupied electron number 
caused by the core-hole potential within the muffin-tin 
sphere at the core-hole site. 
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Figure 3: The d-DOS at the site of 3s-core hole with (a)up 
spin and (b)down spin. 



Up spin core-hole 

Figure 02(a) shows the d-DOS at the site of core hole 
with up spin. By comparing with the d-DOS with no core 
hole, we notice that the attractive core-hole potential is 
strongly working on the 3d electrons with down spin, 
pulling down the d-DOS. A large weight develops around 
the bottom of the band. Even quasi-bound states seem to 
be formed. It looks like an overscreening (the screening 
charge reaches 2.38). On the other hand, the 3d electrons 
with up spin are prevented from coming close to the core 
hole by the exchange effect, thereby the attractive po- 
tential would be less effective. The repulsion between 
3d electrons pushes the 3d states upward, resulting in a 
compensation of the overscreening (the screening charge 
-1.44). As a total, the screening charge is 0.94. This in- 
dicates that the screening is almost complete within the 
muffin-tin sphere on the core-hole site. 



Down spin core-hole 

Figure 02b) shows the d-DOS at the site of the core 
hole with down spin. The attractive potential is now 
strongly working on the up spin states, pulling down the 
d-DOS. Again a large weight develops around the bottom 



of the band. The difference is that the screening charge is 
limited to 0.47 in the up spin 3d states, even after nearly 
all occupied states are pulled down below the Fermi level. 
Therefore, the screening has to be completed by the down 
spin 3d states. Since the attractive potential is less effec- 
tive on the minority spin electrons due to the exchange 
effect, there is no appreciable weight around the bottom 
of the band. The screening is completed only through 
a moderate transfer of weight from unoccupied levels to 
the occupied levels. 



2. XPS spectra 

We construct the excited states by distributing N e 
band electrons on the one-electron states in the presence 
of the core hole. Since the first BZ is quite small, we pick 
up only the T-point and distribute 8 x 27 electrons on 3d, 
4s, 4p bands. The excited state with the lowest energy is 
given by piling one-electron states up to the Fermi level. 
Then, picking up an occupied electron from these states 
and placing it on a some unoccupied state, we create one 
electron-hole pair. We create two electron-hole pairs in 
a similar way. We calculate Eq. (CQ) with these excited 
states by evaluating the Slater determinant in Eq. ((4]) 
with the use of the APW wavefunctions. Since the first 
BZ is quite small, we take account of the states at the T- 
point only. Figured] shows the calculated XPS spectra as 
a function of the binding energy u> q — e in comparison with 
the experiment^ The contributions from zero, one, and 
two electron-hole pairs are separately shown. The calcu- 
lated spectra reproduce well the experimental data. In 
the following, we discuss the origin of the spectra. 



Up spin core-hole 

The excited state with the lowest energy (no electron- 
hole pair), denoted as |/o), has a finite overlap with the 
ground state \g), giving rise to intensities on the thresh- 
old (broadened by a Lorentzian function with FWHM 
2r s = 2.0 eV). The energy of threshold is adjusted to 
coincide with the experimental value. The overlap is 
given by {fo\s^\g) = A T Aj. with A T = 0.980 + i 0.079 
and Ai — —0.055 — i 0.352, where Af and A± represent 
the overlap with respect to up spin electrons and with 
respect to down spin electrons, respectively. In princi- 
ple, such overlap tends to be zero with N e — > oo, ac- 
cording to the Anderson orthogonality theorem.— In such 
infinite systems, energy levels become continuous near 
the Fermi level and thereby infinite numbers of electron- 
hole pairs could be created with small excitation ener- 
gies, leading to the so called Fermi edge singularity in 
the XPS spectra. The finite contribution obtained above 
arises from the discreteness of energy levels and could be 
interpreted as the integral intensity of singular spectra 
near the threshold, according to the model calculations 
to other systemsi&ii Note that \A^\ 2 = 0.966, which is 
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Figure 4: XPS spectra as a function of binding energy for 
the 3s-core hole with (a) up spin and (b) down spin. The 
experimental data are taken from Ref. [la ]. 



close to 1 , in spite of the large modification of the local 
3d-DOS at the core-hole site (see Fig. OJa)). This may 
be understood as follows; since the 3d bands with up 
spin are almost occupied in the ground state, the wave- 
functions of the 3d bands in the presence of core hole 
could be nearly represented by a unitary transform of the 
wavefunctions of the occupied 3d electrons in the ground 
states. Therefore the value of the Slater determinant in 
the up spin channel could not be changed by the pres- 
ence of the core hole. On the other hand, |j4J 2 = 0.127, 
suggesting that the wavefunctions of the occupied lev- 
els are considerably modified from those of the ground 
state in the down spin channel, which modifications are 
brought about by the mixing of the unoccupied states of 
the ground state, because the 3d bands are only partially 
occupied in the ground state. This change reduces the 
value of the Slater determinant in the down spin channel. 

For the excitation of one electron-hole pair, the con- 
tributions from the up spin channel are negligibly small. 
Since the wavefunctions of the unoccupied levels are al- 
most orthogonal to the wavefunctions of the 3d bands 
of the ground state, the Slater determinant correspond- 
ing to the the occupation of such states would be almost 
vanishing. On the other hand, the excitations of one 
electron-hole pair in the down spin channel could give 
rise to a considerable contribution to the spectra; each 
one-electron wavefunction of occupying electrons is mod- 



ified from those of the ground state by the mixing of 
the unoccupied states, resulting in non-vanishing Slater 
determinants. As a result, we obtain the intensities dis- 
tributed with wide range of binding energy, correspond- 
ing to various kind of one electron-hole pair. They consist 
of two peaks; one is located around ui q — e = 92.5 eV and 
another is around 95 eV. The latter peak corresponds to 
the excitation from the quasi-bound states around the 
bottom of the 3d band in the down spin channel. 

For the excitation of two electron-hole pairs, the con- 
tributions from the up spin channel are extremely small 
by the same reason as for the case of one electron-hole 
pair. The contributions arise from the down spin chan- 
nel; they consist of three peaks, one small peak around 
93.2 eV, and large peaks around 96 and 99 eV. The exci- 
tation of three electron-hole pairs gives rise to only minor 
contribution to the spectra. 



Down spin core-hole 

The excited state with the lowest energy (no electron- 
hole pair) has an overlap with the ground state, 
(/oNs) = A r A t with A T = 0.708 + i 0.684, A L = 

0. 758 + i 0.576. Note that |,4 T | 2 = 0.969, which is close to 

1, although the local 3<i-DOS is strongly modified from 
that in the ground state (see Fig.[U[b)). This means that, 
by the same reason as explained in the up spin core-hole, 
the one-electron states in the presence of core hole are 
nearly represented by a unitary transform of occupied 3d 
states of the ground state, and thereby the corresponding 
Slater determinant is close to 1. In the down-spin chan- 

0.906, which is less than \A^\ 2 but close to 
This simply means that the modification of the one- 



nel, |AJ 2 
1 

electron states by the core-hole potential is small, due 
to the exchange potential. Thereby we obtain a large 
spectral weight around the threshold. 

For the excitation of one electron-hole pair, we again 
have extremely small contributions in the up spin channel 
because of the vanishing overlap with the ground state. 
In the down spin channel, since the modification of one- 
electron states by the core-hole potential is small, the 
creation of electron-hole pair makes the corresponding 
state nearly orthogonal to the ground state, and thereby 
small contributions come out, as shown in Fig. [H 



B. ls-core hole 

1. one-electron states 

Figured] shows the d-DOS at the ls-core hole site. It is 
a little modified from that with no core hole (Fig. [2]) . It 
depends little on the spin of the core hole, indicating that 
the exchange effect is rather small. This is anticipated 
from the fact that the Is core is much localized around 
the origin than the 3s core. In addition, no quasi-bound 
state exists around the bottom of the 3d bands. Table 
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Table II: Screening charge with the d symmetry within the 
muffin-tin sphere at the ls-core-hole site. 

An dT An d | An d (total) 
up spin hole 0.35 0.74 1.09 
down spin hole 0.36 0.73 1.09 



0.6 



Ull lists the screening charges caused by the core-hole po- 
tential within the muffin-tin sphere at the core-hole site. 
The screening charges are nearly independent of the spin 
of core hole. 
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Figure 5: The d-DOS at the site of ls-core hole with (a) up 
spin and (b) down spin. 



2. XPS spectra 

Figure [6] shows the XPS spectra for the core hole 
with up spin and down spin, respectively. The spec- 
tra are broadened by a Lorentzian function with FWHM 
2r s = 2.0 eV to take account of the life-time broad- 
ening of the Is core. The lowest energy state (no 
electron-hole pair) in the presence of core hole has a fi- 
nite overlap with the ground state; (fo\si\g) — A-^A^ 
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Figure 6: XPS spectra as a function of binding energy for the 
ls-core hole with (a) up spin and (b) down spin. 



for the core hole with up (down) spin, most intensities 
are concentrated near the threshold. Only small intensi- 
ties arise from the creation of electron-hole pair with no 
appreciate satellite. Unfortunately, we don't know the 
experimental data to be compared with. 



IV. CONCLUDING REMARKS 

We have developed an ab-initio method to calculate the 
inner-core XPS spectra. We have applied the method to 
the 3s-core XPS in a ferromagnetic metal Fe, where the 
spectra are resolved by the spin of photoelectron. We 
have found that the spectral intensity distributes in a 
wide range of binding energy with a clear satellite for the 
core hole with up spin, while the intensity is concentrated 
near the threshold with no satellite peak for the core hole 
with down spin, in good agreement with the experiment. 
The origin of such behavior has been explained in relation 
to the 3c? band modified by the core-hole potential. 

We have considered only the states on the T point in 
the first BZ to which electrons are distributed in the cal- 
culation of the XPS spectra. Since the first BZ is reduced 
to a small size in a system of supercells, this may not 
cause large errors. With increasing the states to which 
electrons are distributed, one may expect that the overlap 
between the lowest energy state in the presence of core 



7 



hole and the ground state would be reduced, and that the 
contributions from electron-pair creation would increase 
near the threshold, leading to an asymmetric peak near 
the threshold as a function of binding energy. Such be- 
havior has been demonstrated in numerical calculations 
on finite size systems.—^ On the other hand, the struc- 
tures with high binding-energy are expected to be only a 
little influenced by such a refined treatment. In any case, 
to be more quantitative, we need to increase fc-points to 
which electrons are distributed as well as to enlarge the 
size of supercells in the calculation of the XPS spectra. 

Acker et al. pointed out that there is poor correlation 
between the satellite position and the magnetic moment 
on the Fe atom in various alloys and compounds, and 
that the satellite splitting is observed even in some Pauli 
paramagnetspi^ We have interpreted the satellite as aris- 
ing from an excitation from quasi-bound states to unoc- 
cupied states in one-electron states. At present it is not 



clear how the position of quasi-bound state is correlated 
to the magnetic moment in various systems and whether 
the satellite appears even in Pauli paramagnets within 
the present approach. We need further studies to clarify 
these points. 

The present method could be extended to analyze XPS 
spectra from other core holes, for example 2p core in Ni, 
and the x-ray absorption spectral These extensions are 
left in future studies. 
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